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We examine the evolution of an Ising ferromagnet endowed with zero-temperature single spin-flip 
dynamics. A large droplet of one phase in the sea of the opposite phase eventually disappears. An 
interesting behavior occurs in the intermediate regime when the droplet is still very large compared 
to the lattice spacing, but already very small compared to the initial size. In this regime the shape of 
the droplet is essentially deterministic (fluctuations are negligible in comparison with characteristic 
size). In two dimensions the shape is also universal, that is, independent on the initial shape. We 
analytically determine the limiting shape of the Ising droplet on the square lattice. When the initial 
state is a semi-infinite stripe of one phase in the sea of the opposite phase, it evolves into a finger 
which translates along its axis. We determine the limiting shape and the velocity of the Ising finger 
on the square lattice. An analog of the Ising finger on the cubic lattice is the translating Ising 
soliton. We show that far away from the tip, the cross-section of the Ising soliton coincides with 
the limiting shape of the two-dimensional Ising droplet and we determine a relation between the 
cross-section area, the distance from the tip, and the velocity of the soliton. 

PACS numbers: 05.50.-fq, 68. 35. Fx, 68.35.Md, 05.70.Np 



I. INTRODUCTION 

At low temperatures, interfaces separating two broken- 
symmetry ordered phases generally shrink and eventually 
disappear. This coarsening process is very complicated as 
it usually involves the evolution of numerous interfaces. 
The general understanding of coarsening (the emergence 
of the coarsening domain mosaic with a single charac- 
teristic scale growing in a universal manner) has been 
steadily improving over the last forty years [D 12], al- 
though many concrete questions remain unanswered [3]. 
Even in the simplest situation when the two-dimensional 
Ising ferromagnet endowed with a non-conservative dy- 
namics is quenched from the temperature above the crit- 
ical to zero temperature, very few analytical results have 
been established. (One such result describes the ultimate 
fate of the system, e.g. the probability to end up in a 
stripe state another [51 gives the statistics of domain 
hulls in the case of the curvature-driven dynamics.) 

Here we shall examine the evolution of a single inter- 
face. We shall always assume that the Ising ferromag- 
net is subjected to zero-temperature non-conservative dy- 
namics and we also consider the two-dimensional setting 
if not stated otherwise. Even in this situation the evolu- 
tion of a single closed interface is not fully understood. 
(Equivalently, the evolution of a simply connected do- 
main of the minority phase surrounded by the sea of the 
majority phase could not be generally 'solved' [3 HI)- 
The detailed evolution of a closed interface is not partic- 
ularly interesting, but the asymptotic is as it is presum- 
ably universal (independent on the details of the initial 
condition). This is one of the goals of this study. 

To examine the evolution we need a precise description 
of the dynamics. Two most popular non-conservative 
zero-temperature dynamics are the single spin-flip dy- 
namics (the chief example of the microscopic dynam- 
ics), and the time-dependent Ginzburg-Landau (TDGL) 
equation which is the prime example of the macroscopic 



dynamics. A zero-temperature spin-flip dynamics for- 
bids energy raising flips. Glauber's version j6j of the sin- 
gle spin-flip dynamics specifically prescribes that energy 
conserving flips occur at a twice smaller rate than en- 
ergy decreasing flips. For the Metropolis algorithm both 
energy conserving and energy decreasing flip rates are as- 
sumed to be equal. The energy-lowering spin-flip events 
do occur, yet their frequency is asymptotically negligible 
and therefore these subtle differences in the single spin- 
flip dynamics are irrelevant for the problems which we 
shall study. In the following, we set the rate of energy 
conserving flips to unity. 

Even at zero temperature the TDGL equation is dif- 
ficult to analyze since mathematically it is a non-linear 
parabolic partial differential equation 'TJ . Fortunately, in 
the interesting situation when the width of the interface 
is much smaller than its radius of curvature the TDGL 
equation reduces to a much simpler Lifshitz-Allen-Cahn 
(LAC) equation 8J which asserts that the normal ve- 
locity of the interface is proportional to the local mean 
curvature. (The terminology is not yet settled: The LAC 
equation is often termed the AUen-Cahn equation; some- 
times the TDGL equation is called the AUen-Cahn equa- 
tion; also in mathematical literature, people usually talk 
about mean-curvature flows [S].) 

The asymptotic evolution of a single closed interface 
in two dimensions is fully understood in the realm of 
the LAC equation. Indeed, the Grayson theorem [10] 
asserts that the interface approaches to a circle, so at 
the flnal stage of shrinking to a point the interface is 
the circle. Thus with respect to the macroscopic non- 
conservative dynamics the limiting shape is the circle. 
What is the limiting shape with respect to the spin-flip 
dynamics? This question was previously investigated by 
Karma and Lobkovsky [TT1 who exploited the self-similar 
behavior on the late stage of evolution and succeeded in 
reducing the problem to an ordinary differential equation 
which they solved numerically. In Sect.|TT]we use an addi- 
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FIG. 1: (Color online) The quarter of the limiting shape of 
the Ising droplet and the quarter of the circle (the limiting 
shape corresponding to the macroscopic curvature-driven dy- 
namics). The Ising droplet encloses the circle. 



tional trick to reformulate the mathematical description 
in terms of the Stefan problem (more precisely, the dif- 
fusion equation with moving boundaries whose position 
is determined in the process of solution). This Stefan 
problem admits an exact self-similar solution which al- 
lows one to analytically determine the limiting shape of 
the Ising droplet. Figure [T] plots both limiting shapes. 

Infinite interfaces exhibit less universal evolution. In- 
deed, the spectrum of qualitatively different initial con- 
ditions is much more broad, and accordingly there are 
infinitely many different limiting shapes. For instance in 
the situation when the minority phase initially occupies 
the wedge (with opening angle smaller than tt) , the LAC 
equation admits a self-similar solution |12j which gives 
the limiting shape (parametrized by the opening angle of 
the wedge). On the square lattice the most natural non- 
trivial 'wedge' is the corner (the positive quadrant). The 
evolution of the corresponding infinite interface can be 
described in great details thanks to the mapping [T31 |T3] 
on the symmetric exclusion process. Using this mapping 
together with results and methods developed in the stud- 
ies of the symmetric exclusion process (see e.g. [IS] and 
references therein) one can determine both the limiting 
shape and fluctuations p2]. 

Another interesting class of initial conditions corre- 
sponds to semi-infinite stripes. In this case, stripes 
quickly approach the limiting shape (known as the finger) 
which propagates with constant velocity along its axis. 
The shape of the finger corresponding to the macroscopic 
dynamics has been rediscovered a few times as it arises 
in numerous applications, e.g. in modeling of the motion 
of grain boundaries in an annealing piece of metal |17j . 



magneto-hydrodynamic models for the solar fiares |18) . 
dendritic crystal growth [19j . boundary renormalization 
group flows [5D], and various other problems in physics 
|21) and mathematics [5^ . For spin-flip dynamics (in this 
case we consider the finger which is parallel to one of the 
axes of the square lattice) the shape of the Ising finger 
was unknown. We compute the shape of the Ising finger 
in Sect, ini 



In Sect. IV we discuss challenges in computing the lim- 
iting shapes in three dimensions and look at a particular 
limiting shape corresponding to a translating Ising soli- 
ton. This is an approximately parabolic object which 
moves at a constant velocity along its axis (coinciding 
with the axis of the lattice). The final Sect. [v| contains 
conclusions. 



II. ISING DROPLET 

Any finite domain of one phase in the sea of the oppo- 
site phase disappears in a finite time that scales as the 
square of the characteristic size; in that sense, we quali- 
tatively understand the shrinking of a finite domain. A 
more detailed quantitative understanding emerges in the 
long time limit since after a proper re-scaling, the inter- 
face admits a deterministic limiting shape as long as the 
droplet remains large. In this section we compute the 
limiting shape of the Ising droplet on the square lattice. 

The spin-fiip dynamics is stochastic and this causes a 
number of subtle differences with the macroscopic dy- 
namics (the deterministic LAC equation). To highlight 
one of the differences we note that the macroscopic dy- 
namics always leads to the decrease of the area of the 
simply connected domain 2). To verify this assertion we 
write the LAC equation in the form Wnormai = —DK, 
where K is the curvature and D is the proportionality 
factor which has the dimension of the diffusion constant. 
The area A{t) — aie&[T){t)] then evolves according to 



dA 
~dt 



dsDK 



-2ttD 



(1) 



The last step in Eq. ([ij is the consequence of the Gauss- 
Bonnet theorem. 

Equation ([I]) helps to understand coarsening dynamics 
of two-dimensional Ising [5 and Potts [22] systems, yet 
its deterministic nature disagrees with the microscopic 
spin-fiip dynamics. For instance, the zero-temperature 
spin-fiip dynamics allows area-raising moves. More pre- 
cisely, any domain D can grow up to its rectangle enve- 
lope (which is defined to be the smallest rectangle aligned 
with the axes of the square lattice that contains D). How- 
ever, on average every domain shrinks. Another subtle 
feature of the stochastic spin-flip dynamics is that a single 
connected domain can evolve into a few disjoint domains 
(see Fig. [2]). Apart from pathological initial conditions, 
e.g. those which contain 'tendrils' of width one (i.e. equal 
to the lattice spacing) or strips of width one, or rare sepa- 
rations of tiny drops (Fig. [2]) , such break ups start to play 
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FIG. 3: An illustrative interface rotated by 7r/4 and the cor- 
responding particle configuration. A spin-flip event is shown 
together with the correspondence hop of the particle in the 
symmetric exclusion process. 



FIG. 2: An illustration of a rare emergence of a short-living 
tiny closed interface during the evolution of the finger. 



a significant role only when the droplet becomes compa- 
rable with the lattice spacing; this late stage of evolution 
is clearly stochastic, but it is not interesting. With all 
these caveats we can talk about single domain, ignore 
fluctuations and focus on the limiting shape as long as 
the domain is very large compared to the lattice spacing. 

A few interesting mathematical papers discussed the 
replacement of the LAC equation with isotropic surface 
tension by an equation with anisotropic surface tension 
which is supposed to be a proper macroscopic equation 
corresponding to the spin-flip dynamics We fol- 

low a similar approach, yet our goal is not to fully jus- 
tify the governing equation (this is essentially achieved 
in previous papers), but rather to solve it analytically. 

The initial shape of the droplet is expected to become 
asymptotically irrelevant. For instance, the initial rect- 
angular envelope x Ly can be significantly different 
from the square, Lx{t = 0) ^ = 0), yet asymptoti- 
cally the rectangular envelop of the droplet is approach- 
ing to the square. Denote by 2L{t) the size of this square 
at time t. For concreteness, let the envelope be the square 
< x,y < 2L. Due to symmetry, we can limit ourselves 
to its quarter, < x,y < L; the boundary of the droplet 
thus goes from {x,y) = (0,L) to {x,y) — {L,0). This 
boundary can be represented by a staircase of kinks that 
can pile-up at the same site, but cannot pass through 
each other. It is more convenient to use a representa- 
tion in terms of an exclusion process, that is a collection 
of particles which undergo a random walk and cannot 
occupy the same site. (Both these representations have 
appeared in the literature, see e.g. pTUl4l [27H29] ). 

The representation in terms of the symmetric exclusion 
process becomes evident after rotating counter-clockwise 
by angle 7r/4 around the origin and projecting the bound- 
ary onto the horizontal line (see Fig.jsj). We put a particle 
on the bond (leave the bond empty) if the corresponding 
bond on the interface goes along co-diagonal (diagonal). 



The particles can occupy the lattice sites, with no more 
than one particle per site, and the particles undergo the 
symmetric exclusion process. 

In the long time 'hydrodynamic' limit we employ a con- 
tinuum description. The particle density n{z, t) satisfies 
a diffusion equation 



dn 



on 



(2) 



on a shrinking interval —L{t) < z < L{t). The boundary 
conditions are 



n(-L(<),i) = 1, n{L{t),t)^Q 



(3) 



It is sufficient to seek a self-similar solution which de- 
pends on z and t only through a combination z/L{t), 



n{z,t)^N{Z), Z^z/L{t) 
Plugging Q into ^ we obtain 

- LLZN' = N" 



(4) 



(5) 



where prime (dot) denotes differentiation with respect to 
Z (time t). We now notice that the width L{t) shrinks 
at a rate that is equal to the ffux of particles: 



L 



dn 
dz 



= -iV'(l) 
--L(t) L ^ ' 



Introducing the parameter 6 defined via 
LL = -2b = N'{1) 

we re-writc ^ as 

N" = 2bZN' 
Integrating Q subject to N'{1) = -26 we obtain 

N' ^-2be''''"'' 

Integrating Eq. ([s]) and using A^(l) = we get 

pi 

N{Z) = 26 / dve'"''-'' 
Jz 



(6) 
(7) 
(8) 

(9) 



The boundary condition A'^(— 1) = 1, or equivalently 
N{0) = 1/2, yields 



1 = 46 / dv e'"' 



(10) 



from which b k, 0.3051025211. In terms of original vari- 
ables, the interface is implicitly given by 



Writing 



X = 



dz n(z, t) 



x-y 



Y = 



L{t)' ^ Lit) 
and using ([9]) we re- write (111 as 



(11) 



(12) 



Y = l-e'"' -'' -2bu I dwe""'-', u = X-Y (13) 



The limiting shape ( 13 



together with the circle (the 
limiting shape corresponding to the LAC equation) are 
plotted on Fig. [T] The Ising droplet encloses the circle 
when both limiting shapes are re-scaled in such a way 
that their rectangular envelops are identical [0, 2] x [0, 2] 
squares. For instance, the diagonal point on the Ising 
droplet lies at 



X* = n = 1 - e 

while for the circle 

X. = K = 1 



1 



V2 



0.262952192433887454 . 



0.292893218813452427... 



Near the axes the Ising droplet is more flat, e.g. 

Y ^bil~ Xf + ... 
when X — > 1, while in the case of the circle 

Y = ^{i-xr+... 

Chayes, Schonmann and Swindle [M] proved that the 
area of the droplet decreases with average rate 4. This 
theorem provides a useful consistency check. Let us first 
compute the area under a quarter of the interface 



dX dY 
Jo Jo 



(14) 



Changing variables, {X, Y) {u = X—Y, Y), and noting 
that the Jacobian is equal to unity, y) = Ij we get 

A^2L^ [ duY{u) (15) 
Jo 
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FIG. 4: Schematic illustration of the evolution of a semi- 
infinite strip (a rectangular finger of the minority minus phase 
surrounded by the majority plus phase). On the right-hand 
side, the flip of the lowest minority spin, the tip spin in this 
example, is an irreversible process that causes the minimum 
height of the finger to advance by one. 



Equation ( 13 ) gives an explicit expression for Y — Y{u) 



Plugging it into ( |15| we compute the integral and find 
that the area of the droplet A = 4L^ — AA is given by 



A = 4L2 



du (1 -I- 6u^)e 



1 



(16) 



Differentiating and using LL ~ —2b we get [3D] 

»i 



A = -166 



du (1 + bu^)e 



- 1 







= -4 (17) 



The last relation is established by using ( 10 ) and 

r-l 



due^" -Hu^ = 1- (46)-i 



(18) 



which is derived from (10 1 through integration by part. 



III. ISING FINGER 

Here we consider the finger geometry (Fig. |4]), e.g. we 
assume that the minority phase initially occupies the 
semi-infinite region y > Q and \x\ < L. The interest- 
ing regime is i ^ L^, where the two corners of the initial 
finger interact and the finger relaxes to a limiting shape 
that eventually recedes at constant velocity. In a refer- 
ence frame moving with the finger, the interface y{x) is 
thus stationary. 

We consider the spin-flip dynamics which is stochas- 
tic, so the area of the finger can occasionally increase, al- 
though on average it decreases. The rectangular envelop 
is now a semi-infinite region y > h and |x| < L where 
h{t) is the current height of the tip. More precisely, the 
tip is formed by all adjacent spins on the lowest height. 
If the tip of the finger contains a single spin (as on the 
right-hand side of Fig. |4|, then when this spin flips, the 
fingertip irreversibly advances by one unit. The finger 
can shed disconnected pieces whenever the tip of the fin- 
ger has the width equal to one and the height greater 
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Ising Finger 
TDGL Finger 





X 



FIG. 5: (Color online) The limiting shape of the Ising finger, 



Eq. ( 22 1, corresponding to the microscopic spin-flip dynamics; 
the limiting shape of the TDGL finger, Eq. (23 1, correspond- 



ing to the macroscopic curvature-driven dynamics. 



than one (see Fig. [2]). Here we consider a very wide fin- 
ger, L ^ 1, and in this situation the above subtleties are 
asymptotically negligible, so we ignore fluctuations and 
shedding events and focus on the limiting shape. (For 
narrow fingers the stochastic effects are important. The 
full description of the evolution of the finger of the least 
possible width 2 is highly non-trivial and unknown.) 

Due to symmetry we can limit ourselves to the region 
< a: < L and y > 0. The governing equation for the 
interface shape y{x,t) is [11] 



yt 



Vx 



(19) 



where Vt = 'g(, Vx = etc. In an upward moving 
reference system in which the surface is stationary 



Vx 



{i + yxy 



(20) 



This equation must have a solution satisfying the bound- 
ary conditions y(0) = and y{L) ~ oo. These conditions 
fix the velocity 



1 

L 



(21) 



The shape of the finger is 



Y 



Hi-x)-x, (x,r) = (|,|) (22) 



This theoretical prediction perfectly agrees with previous 
simulation results (cf. Fig. [s] with Fig. 8 of [T^]) 



It is interesting to compare the Ising finger ( 22 ) with 



the finger that arises if the system evolves according to 
the TDGL equation, that is the mean-curvature evolution 
(Fig. [5]). In the latter framework, the interface y{x) sat- 
isfies the equation y^x = v{l + y'^) in the upward moving 



reference system in which the surface is stationary. Inte- 
grating this equation and imposing the boundary condi- 
tion y ^ 00 when |a;| — > L, one arrives at the following 
TDGL finger 



Y 



In 



ttX 



(23) 



This finger recedes at a constant velocity that is given 
hy V = y"{x = 0) = i:/2L. The expression ([23]) for the 
finger shape has been derived by MuUins in the context of 
modeling the motion of grain boundaries |17| , and it has 
subsequently arisen in numerous applications [T^ [TBl - 
I21j : this finger also appears under the name grim-reaper 
in the mathematics literature and hair-pin j20j in the 
physics literature. 



IV. ISING SOLITON 

In three dimensions, even in the realm of the LAC 
equation (equivalently the mean-curvature evolution) lit- 
tle is known about the limiting shapes. If the initial do- 
main is convex and compact, it approaches to a ball and 
eventually shrinks to a point. Thus in this setting the 
limiting shape is a sphere |31j and the behavior is similar 
to the behavior observed in two dimensions [TU]. How- 
ever, in three dimensions the assumption of convexity is 
important — a compact closed interface which is topolog- 
ically a sphere may become singular before it shrinks to a 
point. The simplest example is provided by a dumbbell 
with sufficiently thin neck. The classification of possi- 
ble limiting shapes (the surfaces of the shrinking LAC 
droplets) and the description of singularities of closed in- 
terfaces undergoing the mean-curvature evolution is an 
active research area, see e.g. Refs. 

In the case of zero-temperature spin-flip dynamics, not 
a single limiting shape is known in three dimensions. In- 
deed, even an equation governing the evolution of the in- 
terface in three dimensions is unknown. In the simplest 
case when when the initial domain is convex, or almost 
convex, the limiting shape should be closed to the sphere. 
Simulations [llj confirm this expectation. 

Consider now possible analogs of the two-dimensional 
finger. Let us first discuss the mean-curvature evolution. 
In three dimensions there are various possibilities, e.g. 
one can start with a semi-infinite bar, or consider rota- 
tionally invariant initial conditions. The former case is 
more similar to the two-dimensional finger, but math- 
ematically the problem is challenging since even in the 
reference frame moving with the finger one must solve a 
non-linear partial differential equation. In the rotation- 
ally invariant setting the problem reduces to an ordinary 
differential equation. The limiting shape is very differ- 
ent from the finger, viz. the width of the cross-section 
diverges with the distance from the tip. This solution 
is usually termed the (rotationally invariant) translating 
soliton and it has been investigated in Refs. i34-364. Us- 
ing the cylindrical coordinates and representing the inter- 



6 



face in the form r = r(z, t) we reduce the mean-curvature 
evohition equation to 



role of R, cf. (25). Thus the shape of the Ising sohton 



rt 



l + ri 



(24) 



This equation admits a family of solutions parametrized 
by the soliton speed u > 0. For the soliton moving with 
constant speed v along the z— direction we write 

r(z, t) = v-^R{Z), z-vt = y-^Z (25) 

and recast Eq. ( |24[ ) into an ordinary differential equation 



R' 



R" 



1 + iR 



■A2 



1 

R 



(26) 



where prime denotes the derivative with respect to Z. 
Equation (26) cannot be solved in quadratures, but one 



can readily extract asymptotic behaviors near the tip of 
the soliton (without loss of generality we set the position 
of the tip to Z = 0) and far away from the tip. Near the 
tip, Z +0, one finds that R^ has a regular asymptotic 
expansion in powers of Z: 



R^ ^4:Z- \Z'^ - ^Z^ 



(27) 



Far away from the tip, Z — )■ oo, the expansion is more 
cumbersome as it contains logarithms; 



R' 



2Z + \nZ 



(28) 



Overall, the interface is approximately parabolic. 

Consider now an Ising soliton on the cubic lattice 
which moves with a constant velocity along z axis. This 
Ising soliton is an analog of the rotationally invariant soli- 
ton (26)-(28). We cannot provide the detailed descrip- 



tion of the Ising soliton as we don't even know the govern- 
ing evolution equation. Far away from the tip, however, 
the problem is essentially two-dimensional. (The same 
is valid in the case of the macroscopic dynamics — far 
away from the tip the LAC equation ( 24 ) reduces to the 
two-dimensional LAC equation rt = —1/r.) Therefore 
the cross-section of the Ising soliton must be identical 
to the Ising droplet which was studied in Sect. |ll| Let 
[—L,L] X [—L,L], where L — L{z,t), be the rectangular 
envelop of the Ising soliton. The soliton moves without 
changing its shape, so L{z,t) — L{Z) with Z defined in 
( [25] ). Plugging L{Z) into ^ we get 



dL 
'dZ 



-v'^L^ = -2b 



from which we find the leading behavior far away from 
the tip: 



(29) 



This behavior is similar to (28) describing the rotation- 



ally invariant soliton once we notice that vL plays the 



is asymptotically parabolic, e.g. the area of the cross- 
section scales linearly with the distance from the tip in 
the Z — >■ oo limit; the cross-section is not a disk even far 
away from the tip, it is actually asymptotically identical 
to the Ising droplet on the square lattice. 

To determine the entire shape of the Ising soliton one 
must know an equation governing the evolution of the in- 
terface in three dimensions and even then one would rely 
on numerical integration since even in the reference frame 
moving with the soliton a governing equation would turn 
into a non-linear partial differential equation. One pos- 
sible direction is to guess a governing equation and then 
try to determine the behavior near the tip where one 
can employ asymptotic methods. An attempt of such an 
analysis is given in Appendix [Xj 



V. CONCLUSIONS 

We have performed an analytical study of the limit- 
ing shapes arising in the context of the zero-temperature 
single spin-fiip dynamics. We have determined the limit- 
ing shape of the shrinking Ising droplet and the limiting 
shape of the Ising finger which moves along its axis (coin- 
ciding with one of the two axes of the square lattice) . We 
have investigated the limiting shapes only on the square 
lattice; an analytical computation of the limiting shapes 
on the hexagonal lattice appears feasible and it would be 
interesting to perform such a calculation. 

There is much more room for diverse limiting shapes 
on the cubic lattice. Nothing is known about these lim- 
iting shapes. The chief reason is the lack of the equation 
describing the evolution of interfaces on the cubic lat- 
tice. There are also intrinsic mathematical difficulties in 
analyzing such would-be governing equations. For one 
particular limiting shape, the translating Ising soliton 
which moves along its axis (coinciding with one of the 
three axes of the cubic lattice), one can circumvent the 
aforementioned challenges and extract partial analytic in- 
formation. We have shown that asymptotically (that is, 
far away from the tip) the surface of the translating Ising 
soliton is parabolic, e.g. the cross-section area grows lin- 
early with the distance from the tip. Furthermore, the 
cross-section of the Ising soliton coincides with the lim- 
iting shape of the Ising droplet on the square lattice. 
Making a plausible guess about the governing evolution 
equation, we have deduced the shape of the translating 
Ising soliton near its tip (Appendix [A|). 
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Appendix A: Tip of the Translating Ising Soliton 



An additional transformation 



We need to know the equation governing the evolution 
of the interface on the cubic lattice. In a related (but sim- 
pler) problem of spin-flip dynamics in a magnetic field, 
the governing evolution equation for the interface has 
been guessed in recent study [37 . This guess has been 
guided by symmetry considerations (the governing equa- 
tion for z(x, y; t) should be invariant under the change 
of any coordinate pair) and by the requirement that it 
should reduce to the two-dimensional evolution equation 
in appropriate settings. Two functionally independent 
simple equations have been found [37 , as well as families 
of equations built from the independent solutions. Only 
the two independent equations seemed sufficiently simple 
and plausible, one of the two simple equations was found 
to be in excellent agreement with simulations. A similar 
program is possible in the present case, and one gets the 
analog of equation which is expected to be exact is 



(1 



zt = 



(1 



^(1 



~~2 



*^xy 

Zx Zy 



^yy 



(Al) 



This equation is a generalization of Eq. ( 19 ) governing 



the evolution of the interfaces on the square lattice and 
it is an analog of the equation describing the evolution of 
the interfaces on the cubic lattice in the presence of the 
magnetic field |37|. 



Equation ( Al ) applies to the region x > 0,y > where 
Zx > 0,Zy > 0. Generally the first derivatives should be 
written as \zx\ and \zy \ and then the equation will be ap- 
plicable everywhere. Due to symmetry we can limit our- 
selves to the region a; > 0, y > 0, so the form ( Al ) suffices. 



Equation ( Al ) is harder to judge than equations analyzed 



in [37] as we cannot solve (Al) analytically. Postponing 
a careful consideration of the agreement between the pre- 
dictions of ( |A1[ ) and numerical results for the future, let's 
just examine the behavior near the tip. First, we notice 
that for the Ising soliton which moves with constant ve- 
locity V along z axis Eq. (Al) reduces to 



2 / 



X — V ^X, y — y z — vt — V 



(A3) 



allows us to recast Eq. (A2) to f— independent form 



(1 



Zy I 



z 



XX 



(1 



Zx +Zy 



72 



Z^XY 

ZxZy 



Z- 



YY 



(A4) 



Near the tip, that is, when < X < 1 and < F < 1, 
we seek the solution as an expansion in the form analo- 



gous to (27), namely 



(A5) 

Plugging (|A5|) into (|A4|) and keeping only the leading 



Z = C{X^ + 2\XY + 



terms we find the consistency when an 'asymmetry' pa- 
rameter A is the root of the cubic equation A'^ — 3A-I-2 = 
and the amplitude is given by C = 5 (1 + A)'^. Since 
A''' - 3A + 2 = (A - 1)2 (A + 2) there are two roots, A = 1 
and A = — 2. The second root is inappropriate since the 
resulting quadratic form A"^ — AXY + Y^ vanishes in the 
X > Q^Y > Q. Hence A = 1 and C = 2, so near the tip 



Z = 2{\X\ + \Y\f 



(A6) 



where we have written the solution in the form which is 
vahd for all X,Y satisfying |A:| < 1 and |F| < 1. More 
precisely, ( A6 ) is valid when 



v<|X|<l, w<|r|<l 



(A7) 



Indeed, the continuum description underlying the usage 
of evolution equations like ( |Aip applies only when the 
distance from the tip far exceeds the lattice spacing. In 
other words, ^ 1 and \y\ ^ 1, and this in conjunction 
with (A3 1 gives the lower bounds in (A7). Not surpris- 



ingly, the velocity must be very small. The same applies, 
of course, to the two-dimensional finger, viz. the con- 



tinuum description leading to the limiting shape ( 22 ) is 



(A2) valid when the corresponding velocity (21 1 is very small 
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